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ABSTRACT 

In recent years there has been some progress towards detecting solar-like oscillations 
in stars. The goal of this challenging project is to analyse frequency spectra similar 
to that observed for the Sun in integrated light. In this context it is important to 
investigate what can be learned about the structure and evolution of the stars from 
such future observations. Here we concentrate on the structure of the upper layers, 
as reflected in the phase function. We show that it is possible to obtain this function 
from low-degree p modes, at least for stars on the main sequence. We analyse its 
dependence on several uncertainties in the structure of the uppermost layers. We also 
investigate a filtered phase function, which has properties that depend on the layers 
around the second helium ionization zone. 
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1 INTRODUCTION 

Since the successful development of helioseismology in the 
last few decades, the interest in extending the work to other 
stars has become clear. However, the small amplitudes ex¬ 
pected in the power spectra of stars with solar-like oscilla¬ 
tions makes this a very difficult task and to date no unam¬ 
biguous detection has been made. The main limitation is the 
atmospheric noise. 

On the other hand, these stellar spectra are expected 
to have regular patterns that can be detected even with¬ 
out a full determination of the p-mode frequencies. Some 
efforts have been made in this direction (Brown et al. 
1991; Pottasch, Butcher & van Hoesel 1992; Kjeldsen et 
al. 1995). The most interesting parameters that can be ob¬ 
tained in this limit are the so-called small and large separa¬ 
tions, as demonstrated by several theoretical analyses (e.g. 
Christensen-Dalsgaard 1984, 1988, 1993; Ulrich 1986; Gough 
1987; Gough & Novotny 1993). 

However, the full power of asteroseismology for solar¬ 
like stars would, as in the case of helioseismology, require 
extensive determination of individual p-mode frequencies. 
In fact, once the observed amplitudes are above the noise 
level, a few weeks of observations will be enough to obtain 
accurate frequency measurements. In the near future there 
will be space missions devoted to asteroseismology (Baglin 
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1991; Catala et al. 1995) which, we hope, will give the first 
firm results on this issue. With these future observations 
in mind, in this paper we shall assume that a full set of 
low-degree (l < 2) p-mode frequencies are available in a 
given frequency range. We do not consider modes of higher 
degree because they cannot be detected with simple photo¬ 
metric observations. Also, we suppose that the radial order 
n and the degree l of the modes are known, and, of course, 
that the frequencies have been corrected for their rotational 
splittings. The determination of n and l requires compar¬ 
ison with the models but for solar-like stars on the main 
sequence, such as we shall consider here, it is plausible to 
assume that this is possible without any uncertainty. 

As in the solar case, a direct comparison of theoreti¬ 
cal and observed frequencies would be inconclusive. How¬ 
ever, stellar acoustic oscillations satisfy a simple asymp¬ 
totic relation which allows the separation of the contri¬ 
bution of the upper layers, such as the convective enve¬ 
lope, from that of the deep interior. This relation has 
been used extensively in the solar case, for instance to ob¬ 
tain the sound speed (e.g. Christensen-Dalsgaard, Gough & 
Thompson 1989) and to analyse the upper layers (e.g. Perez 
Hernandez & Christensen-Dalsgaard 1994b). Since for the 
distant stars we expect to detect only low-degree modes, we 
re-analyse the relation in this context. A similar investiga¬ 
tion was carried out recently by Lopes et al. (1997), to inves¬ 
tigate what would be the diagnostic potential of observing 
solar-like oscillations in the star (3 Virginis. 

Analysis of observations of low-degree solar p-modes is 
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illustrative of what can be inferred for other stars. Thus we 
first present an analysis of the observations by Lazrek et al. 
(1997), obtained with the GOLF instrument on the SOHO 
spacecraft. Having demonstrated that it is possible to carry 
out the asymptotic analysis for stars, we concentrate our 
work on the upper layers as probed by the so-called phase 
function. We investigate the information that can be ob¬ 
tained from this function, taking into account the plausible 
accuracy of the frequency determinations. We pay particular 
attention to the second helium ionization zone. 


2 THE SUN AS A STAR 


2.1 The phase function for low-degree modes 


Neglecting the perturbations in the gravitational potential 
and considering the waves locally as piano-parallel under 
constant gravity, th e frequencies of p mode s satisfy the 
asymptotic relation (Deubner fe Gough 198^) 


7r (n + e) 
CJnl 



(i) 


Here ui is the angular frequency, c the adiabatic sound speed, 
N the buoyancy frequency, ui c the cut-off frequency and Si 
the Lamb acoustic frequency. The integral is over radius r 
with the limits ri and r 2 defined by the vanishing of the 
bracket in the integral. In principle, the asymptotic theory 
gives e = —1/2 but since close to the surface the asymptotic 
conditions are not satisfied, e is assumed to be a function of 
frequency, yet to be determined - for low-degree modes its 
dependence on l is negligible. 

It is interesting to note that with the sole assumption 
that the perturbations in the gravitational potential can be 
neglected, equation (|lj) still holds in the asymptotic limit, 
provided ui c a nd N are replaced by more general functions 
(Gough 1996). 

For the Sun, and many other stars, equation (|IJ) can be 
further approximated to yield a very useful relation. If there 
is a point ro - usually in the convective envelope - such that 
|lV 2 (r)| <c u> 2 and u> 2 (r) <C oj 2 for r < r o and Sf <g ui 2 
for r > ro, then it can be shown that equation (ij) can 
be approximated by (e.g. Christensen-Dalsgaard & Perez 
Hernandez 1992) 


— ^ F(w/L) - GO , 

Mnl 


where 
F(w) = 
and 

GO = 




c 2 \ 1/2 dr 

o o I 


ira(u) 

u 


( 2 ) 

(3) 

(4) 


Here L = l + 1/2, l being the degree of the mode, and 
w = w/L. The integral in equation (^) is from the lower 
turning point r t (defined by the vanishing of the bracket 
in the integral) to the surface radius R. Finally, the phase 
function a(u/) depends on conditions near the stellar surface. 

This phase function can be computed from the structure 
of the upper layers of a stellar model, by fitting numerically 
computed eigenfunctions to an asymptotic approximation at 



Figure 1. The continuous line is Gas A') for Model A. The dashed 
line is G(uj ) for the same model and the frequency set indicated 
in the text. The dot-dashed line is G cow ©) for the same model 
and mode-frequency set but in the Cowling approximation. The 
dotted line is G + d/u 2 , where d is a constant. The functions have 
been shifted by constants to match at the lowest frequency. 


a point where the latter is valid (see Christensen-Dalsgaard 
& Perez Hernandez 1992). We denote the phase function 
obtained in this way a as (w) and the corresponding related 
function as defined in equation (Q) is denoted G as ( u>). Since, 
in this case, the wave equations are integrated from the sur¬ 
face down to a point in the envelope, by construction a as 
and Gas depend only on the upper layers. The continuous 
line in Fig. I] corresponds to GasO for a solar model with¬ 
out diffusion, (Christensen-Dalsgaard, Proffitt & Thompson 
1993), in the following referred to as Model A. 

It is interesting to note that the well-known Tassoul 
equation (Tassoul 1980), which is an approximation valid 
for low-degree modes, 

(5) 

Vnl 

where u n i is the cyclic frequency and 6, 5, A and A v are 
constants, is a particular case of equation (jiij). In fact, it is 
straightforward to show that in this approximation 

2tt 2 AAz/ 


v n i~ (n + | + o'j Av — ( AL 2 — 5) 


F(w) ~ 


and 


2A v 2w 


+ 


( 6 ) 


. 7T@ 2it 2 5Au 

GO ~ — +-j— 

UJ UI 2 


(7) 


A detailed derivation of equation (|^) from equation (^) was 
given by Vorontsov (1991) (see also Gough 1986). 

If the radial order n and the degree l are assumed to 
be known, F(w) and G( ui) can be estimated from p-mode 
frequencies by fitting them to equation (^|). Hence these 
functions are potentially observable quantities for solar-like 
stars. This is the advantage of using equation (j|2j) rather than 
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the formally more precise equation (jllj). In the present work 
we have made a least-squares fit, expanding F(w) and G( ui) 
in terms of Legendre polynomials in w and ur. 

k m 

— ~a 0 +V'a i / i (a;) +V'feiSi(y) , (8) 

Unl t -— 1 

i= 1 i=l 

where fi and gi are Legendre polynomials of order i with 
coefficients a; and 6;. The variables x and y are linearly 
related to w and uj, respectively, and are defined in the in¬ 
terval [—1,1]. Note that from such a fit F(w) and G(o>) can 
be obtained only to within a constant. 

Throughout the paper we use the same mode data set, 
consisting of modes with l = 0,1, 2 and 12 < n < 27. The 
interval in n has been chosen such as to reject modes with 
very low amplitudes, which are not expected to be observed 
in solar-like stars. We note that rather similar results are ob¬ 
tained if modes with 1 = 3 are added. The fit was carried out 
through y 2 minimization, using frequency errors obtained in 
early observations with the GOLF instrument on the SoHO 
spacecraft (Lazrek et al. 1997). In equation (^j) we have used 
Legendre polynomials of order 14 for each function. We note 
in passing that, as suggested by equation (jtij), the expansion 
of F(w) is most reasonably carried out in terms of w _1 . In 
this case a value fc ~ 5 in the expansion of F(w) is sufficient 
to achieve a good fit. But in any case this does not affect to 
the determination of G(u>) with which we are mainly con¬ 
cerned in this work. 

This kind of fit for G(u>) is rather similar to that used 
in inversion techniques (e.g. Dziembowski, Pamyatnykh & 
Sienkiewicz 1990) though in that case the fit is usually done 
to a function corresponding to differences between models 
rather than to G(lo) directly. On the other hand, Vorontsov, 
Baturin & Pamyatnykh (1992) applied a very similar fit to 
solar p modes, but considering modes with higher degrees. 

For the Sun, equation (^j) is a good approximation for 
modes of moderate degree, in particular for those modes 
with inner turning points between the base of the convec¬ 
tion zone and the second helium ionization zone. When this 
kind of mode set is used, the asymptotic function G as (aj) 
and the numerically fitted G(lj) agree very well. Therefore, 
the observational G( u>) can be used as a test of the upper 
layers. However, when only low-degree modes are included, 
such as is the case for other stars, the approximations lead¬ 
ing to equation (^|) are not so accurate and hence larger 
differences between G as and G are expected. This can be 
seen in Fig. [j], where the dashed line corresponds to G( u>) 
for Model A (since G is obtained to within a constant, in 
Fig. |l] and the following figures we choose the constant such 
that all the functions agree at a given frequency). Indeed, 
there are significant differences between G{uj) and G as (u;), 
the main difference coming from neglecting the perturba¬ 
tions in the gravitational potential in equation (Q). This can 
easily be checked by determining the corresponding func¬ 
tion Gcow(w) by fitting to p-mode frequencies computed 
in the Cowling approximation. This is shown by the dot- 
dashed line in Fig. []J It is indeed much closer to G as . Thus, 
the function G obtained from the observations is a sum of 
G aa ~ Gcow that depends on the upper layers and a func¬ 
tion that depends on the perturbation in the gravitational 
potential throughout the solar interior. 


For completeness we also note that the function F(w) 
obtained from a fit to equation (j^) does not agree with the 
asymptotic expression (LI), except if the frequencies in the 
Cowling approximation are used. 


2.2 Perturbations on the gravitational potential 

In principle, equation (JSj) can be improved by including the 
perturbation in the gravitational potential. To first order, 
the result is an additional term with a specific dependence 
on ui and L, such that equation (|2) is replaced by 




(9) 


(Vorontsov 1991). It is straightforward to show that if F ( w ), 


G(uj) and P$(w) are solutions of a fit to equation (jOj) then so 
are F(w) +C, G(ui) +C+ D/lo 2 and P${w) + D for any con¬ 
stants G and D. In particular, it is not possible to separate 
the contributions of the upper layers, as given by G as , from 
the zero-order term of P$ since in any fit to observational 
data both would be included in G(w). 

We consider again the function G obtained from a fit to 
equation (|2i). In Fig. [ljthe dotted line represents G{uj)+d/uj 2 
with a suitable (positive) value for the constant d (as well 
as including the constant shift mentioned above). It is clear 
that G ~ G as — d/u> 2 . Had we done a fit to equation ([]), 
we would have obtained a similar discrepancy due to the 
undetermined constant D. Thus, in what concerns the de¬ 
termination of G, a fit to equation (R) is equivalent to a fit 
to equation (g). 

We shall now argue that, when considering a pair of 
models (or a model and the observations), the differences in 
P# are very small and hence 5G ~ <5G as . To a first approx¬ 
imation, the function Pj> is given by (e.g. Vorontsov 1991; 
Christensen-Dalsgaard 1996) 






( 10 ) 


where p is the density. Since we are dealing with low-degree 
modes, Pj> can be approximated by a Taylor expansion in 
w = L/oj around the centre. From equation (|I^), we have 

pR j /*1 ~ 

P*(w = 0) = 2 ttG / p — oc p 1 ^ 2 / f;&x , (11) 

Jo c Jo c 

where p is the mean density, x = r/R and p, c are dimension¬ 
less functions. Also, by setting 2 = r 2 /c 2 as the independent 
variable and integrating equation (llOj) by parts twice, it fol¬ 
lows that 


dPj, 

dui 


= 0 . 


w=0 


( 12 ) 


Thus, the variation of P<s> is of second order in w, and in the 
w interval spanned by our data set P$ can be estimated by 
its value at the centre. 

From equation ( ^TJ) , it follows that, for models with 
the same mean density and similar structure in the inner 
region, the functions Pj> are similar; hence, when consider¬ 
ing differences between pair of such models, this term can 
be neglected. In fact, numerical tests show that for realis¬ 
tic solar models the differences between G( u>) are very close 
to the corresponding differences in G as (cu). Hence, although 
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Figure 2. The function G(uj) for the observations (continuous 
line with error bars), Model A (dashed line), Model B (dotted 
line) and Model C (dot-dashed line). The functions have been 
shifted by constants to match at the lowest frequency. 


For distant stars with poorly known parameters, the 
uncertainty in P$( 0)/u> 2 can eventually become comparable 
with the differences in G as . In this case, it may be better to 
work with a phase function invariant under the transforma¬ 
tion G' —* G+C+D/u) 2 or, equivalently, a' —* a+CuJ+D/ uj. 
It is easy to show that the phase function 


7 


da 

= a — uj— - u> 

d oo 


2 


d 2 q 

duF 


(13) 


has this property and hence depends only on the upper lay¬ 
ers. This is a generalization of 


(3 = a — ui 


da 
d ui 


(14) 


which has been used in earlier investigations (e.g. Brodsky & 
Vorontsov 1988; Christensen-Dalsgaard & Perez-Hernandez 
1991; Lopes et al. 1997), and which is invariant only under 
the transformation a —> a + Cu>. 

We note that the function 7 introduces a filter which 
suppresses partially the contribution of the uppermost lay¬ 
ers. Thus, if a good estimate of the uncertainty in P# as 
given by equation (0 is available, it is generally better to 
work with a(u>) or G(ui). On the other hand, if we are inter¬ 
ested in deeper layers (such as the second helium ionization 
zone), then the filtering procedure described below may be 
preferable, in isolating more effectively the contribution from 
just these layers. 


the function G(u>) depends also on properties of the inte¬ 
rior, differences between such functions depend mainly on 
the differences between the uppermost layers. For distant 
stars the mean density is generally not known with suffi¬ 
cient accuracy; however, in this case it is likely that similar 
cancellation can be achieved by working in terms of fre¬ 
quencies measured in units of the large frequency separation 
Au n i = Vni — Vn-11 (cf. equation]^) which scales approxi¬ 
mately as p -1 / 2 . 

In order to illustrate the dependence of G(w) on stellar 
structure, in Fig. ^ we show this function for three solar mod¬ 
els and the mode set indicated above. The dashed line cor¬ 
responds to the previously mentioned Model A, the dotted 
line to a model that includes helium settling (Christensen- 
Dalsgaard, Proffitt & Thompson 1993, hereafter Model B) 
and the dot-dashed line to a solar model that in addition 
has an artificial increase in the surface opacities in order 
to get a better agreement with the observations (hereafter 
Model C). As a result of helium settling, the envelope he¬ 
lium abundance Y env ~ 0.25 of Models B and C is lower 
than the value ~ 0.28 for Model A. On the other hand, 
Model C differs substantially from the other two in the up¬ 
permost layers. As can be seen, the function G(w) for this 
model is also significantly different. In Fig. S we also show 


G(uj) computed with the observed frequencies (Lazrek et 
al. 199^), with errors estimated by Monte-Carlo simulation. 
As expected, Model C is closest to the observational data. 
(Note that, since the superficial layers of the solar model 
suffers from other uncertainties, such as non-adiabatic and 
convective effects, it cannot on this basis be concluded that 
Model C is better than the other two.) This result is similar 
to that obtained by using intermediate-degree modes. 


2.3 The filtered phase function 

Although G depends mainly on the uppermost layers, it also 
has a small contribution from deeper layers. In fact, it is pos¬ 
sible to obtain from G(w) a function which predominantly 
reflects the properties of the layers around the second helium 
ionization zone. 

Perez Hernandez & Christensen-Dalsgaard (1994a) 
showed how a similar function can be extracted through a 
filter. Rather than equation (J2J) , they considered an equiv¬ 
alent expression for differences, assumed to be small, be¬ 
tween pairs of models. With this assumption the differences 
between phase functions are line arly related to differences 
between the equilibrium models (Christensen-Dalsgaard & 
Perez Hernandez 1992) by 


n 2 = 


n8a 


K c (r)—(r) + K^(r)—(r) 

C LO a 


dr 


(15) 


where u> a is the Lamb acoustical cut-off frequency and the 
kernels K c and /© a are known functions. Due to the be¬ 
haviour of the kernels in the solar interior, a localized per¬ 
turbation has a frequency dependence of the form (Perez 
Hernandez & Christensen-Dalsgaard, 1994a) 


H 2 oc cos 2 [ut — (a + l/ 4 ) 7 r] , 
where t is the acoustical depth, 


idr . 

c 


(16) 


(17) 


Hence shorter-period components in 7i 2 correspond to per¬ 
turbations in deeper layers. Thus by filtering the smooth 
components of TL 2 it is possible to obtain a function Tfl 
that depends on deeper layers. 
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Figure 3. Response function H(r ) as function of acoustical depth 
t and radius (on top). The solid line is for equation (^) and the 
dashed line for harmonic functions, both for fits with p = 6. 


In the present case, G(u) can be thought of as the differ¬ 
ence between the actual model (or the Sun) and a smoothed 
model that does not have this localized feature. Relation 
© is then still valid and we can obtain a filtered G f in a 
similar way to Tif. 

The filter used by Perez Hernandez & Christensen- 
Dalsgaard (1994a) was a recursive filter especially suitable 
for H 2 because of the flatness of its smooth component at 
low frequencies (v < 1.4 mHz for the Sun). In the stellar case 
modes at such low frequencies are unlikely to be observed 
because of their expected low surface amplitudes. Further¬ 
more, for the function G, upon which much of our analysis is 
based, the smooth component does not have this behaviour. 
However, below we describe a different filter which has sim¬ 
ilar properties and can be extracted directly from the least- 
squares fit used to compute G(lj). 

Since, in the present work, we obtain G as a Legen¬ 
dre polynomial expansion, a simple way of separating the 
smooth and oscillatory frequency patterns is to consider the 
low- and high-order polynomial coefficients, respectively. A 
very similar type of filtering was introduced by Vorontsov et 
al. (1992). Specifically, we define 

m 

G f (w) = ^ bigi(u>) . (18) 

i=P +1 

A suitable value of p will give the desired oscillatory func¬ 
tion. Similarly, the asymptotic function G as can be fitted to 
Legendre polynomials, from which G aa may be defined. 

In order to choose the value of p in equation (|I^), we 
have computed the response to the filtering of signals of the 
form ©■ Specifically, we define the response function H (r) 
by 


H(r) 


E ,yi(ui,r) 


(19) 


where x(oj,t) is the input signal and y(u>,r) is the output 
function, calculated as in equation (|l^) by expanding x(u, r) 
in Legendre polynomials of order m and setting to zero the 
first p terms; the summation is over the frequency range 
considered. Here we have used the observed frequencies of 
Lazrek et al. (1997). 

We consider filtering corresponding to p = 6. The 
dashed line in Fig. shows the response to simple harmonic 
functions Xi = cos(2 w;t). However, since in equation ( |l6| ) 
a is a function of frequency, the contribution to G(lj) of a 
given layer is not exactly a harmonic function. To investigate 
the effect of this, in equation © we take as input 

Xi = cos2[wir — (a as (a)i) + l/4)7r] , (20) 


with a aa computed for a given solar model. The solid line in 
Fig. I corresponds to H(r ) calculated in this way, with a aa 
obtained from Model A. This term clearly has a significant 
effect, corresponding approximately to a shift in t. From the 
figure, it follows that the filter considered suppresses most 
of the signal from the uppermost layers (r > 0.99717) while 
passing the signal from r ~ 0.9877, where the second helium 
ionization zone is located. Hence in the following we shall 
take p = 6 in equation (fift), at least for the Sun. 

In Fig. |f] we show G^(w) and G f (ui) for a solar model. 
The dashed line corresponds to Gl s , the dot-dashed line to 
Gcow computed in the Cowling approximation and the solid 
line to G f using the full equations. Since G f ~ G aow , G f 
does not depend on the perturbations in the gravitational 
potential. Furthermore, since G f ~ G aa , the function G f de¬ 
pends on the layers between ro (a point in the adiabatically 
stratified convection zone) and an upper limit given by the 
response function, as in Fig. []. 

In Fig. P we show G f (cu) for the same solar models as 
in Fig. The function G f is very similar for Models B and 
C, showing that the filtering has suppressed the substan¬ 
tial difference in G between Models B and C (cf. Fig. ^), 
caused by the differences in the structure of the superficial 
layers. On the other hand G f for Model A is clearly differ¬ 
ent, reflecting the differences in the equilibrium structure 
in the second helium ionization zone. As shown by Perez 
Hernandez and Christensen-Dalsgaard (1994b) this kind of 
filtered phase function depends on three main properties: 
the equation of state, the helium abundance y env in the en¬ 
velope and the specific entropy s in the adiabatically strat¬ 
ified part of the convection zone, whose value controls the 
depth of the convection zone. Model A has the same equa¬ 
tion of state than Models B and C but y» nv is larger and the 
depth of the convection z one smaller. In Fig. we also show 
G f for the observations ( Lazrek et al. 1997| j, with errors 
determined from a Monte-Carlo simulation. It can be seen 
that the models with diffusion (and hence Y env ~ 0.25 and 
a depth of the convection zone closer to that inferred from 
sound speed inversions) agree better with the observations. 
This is in accordance with the the analysis with moderate- 
degree modes, as shown in Perez Hernandez & Christensen- 
Dalsgaard (1994b). Although in the present case the errors 
are larger, our purpose is to show that by using low-degree 
modes alone it is possible to carry out a similar analysis for 
the distant stars. 
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Figure 4. The dashed line is G^ s (ca) for Model A. The solid line 
is G f (aj) for the same model and the dot-dashed line the same 
function but using p-mode frequencies in the Cowling approxi¬ 
mation. 



Figure 5. The continuous line with error bars is G f for the ob¬ 
servations, the dashed line is G f for Model A, the dotted line is 
for Model B, and the dot-dashed line is for Model C. 

3 PHASE FUNCTION FOR STARS 
3.1 Global properties 

In this section we shall consider the phase functions G(ui) 
and G f (u;) for models of stars on the main sequence. The 
models considered are summarized in Table 1. There are 


Table 1 . Stellar models considered. X c is the central hydrogen 
abundance 


model 

M/Mq 

W 

L/L 0 

T e ff 

i 

0.85 

0.6928 

0.329 

5038 

2 

0.85 

0.0003 

0.891 

5466 

3 

1.00 

0.6928 

0.718 

5628 

4 

1.00 

0.1004 

1.238 

5813 

5 

1.30 

0.6928 

2.572 

6509 

6 

1.30 

0.0516 

3.695 

6142 

7 

1.70 

0.6928 

8.512 

8071 

8 

1.70 

0.0479 

10.994 

6619 


two models for each mass, one at the zero-age main sequence 
and the other near the end of hydrogen burning. We shall 
consider the same mode set than for the Sun, that is l = 0,1 
and 2 and 12 < n < 27. To obtain estimates of the errors in 
these functions we assume normally distributed frequency 
errors with a = 0.5 pHz. 

Fig. 1 show G as - dashed line - and G(w) with the 
corresponding errors - continuous line - for all the models 
considered. As for the solar case, most of the differences 
between G as and G can be represented as a function of the 
form d/u > 2 , as shown in the figures (the dot-dashed lines 
correspond to G + d/oj 2 ). We have also computed G(u>) in 
the Cowling approximation. Although for clarity we do not 
show it here, we note that in most of the cases G as — G co w, 
and hence the same comments apply as for the Sun. 

However, for the models at the end of hydrogen burn¬ 
ing, and in particular for the most massive ones, the agree¬ 
ment between G as and G co w is worse, hence indicating that 
other simplifications leading to equation (^|) are not so good 
for these models. In fact, as the stars evolve, the buoyancy 
frequency increases in the core. This effect is more signif¬ 
icant for models with shrinking convective cores, that is, 
for relatively massive stars. As a result, the term in N 2 in 
equation (1) must be taken into account, and hence the sim¬ 
ple asymptotic expression in equation (2) is no longer ade¬ 
quate. Furthermore, the buoyancy frequency introduces ad¬ 
ditional effects in the correction for the perturbation to the 
gravitational potential. [For more complete discussions, see 
for example Vorontsov (1991); Gough (1993); Roxburgh & 
Vorontsov (1994).] Indeed, as shown in Fig. the differences 
between G and G as for the evolved models are not exactly of 
the form d/ui 2 . Thus analysis of data with very small errors 
requires an asymptotic description which is more accurate 
than equation (2). 

To compute the filtered function G f we first look for 
the value of p in equation d) . The goal of the analysis is to 
isolate the signature of the second helium ionization, as was 
done in the solar case in Fig. |E|. Thus the optimal value of p 
must reflect the acoustical depth of the ionization region, as 
indicated by equation (^b|). The location of the second he¬ 
lium ionization varies in temperature from about 1.2 x 10 s K 
in Model 1 to 4.7 x 10 4 K in Model 7; this variation, com¬ 
bined with the change in the stellar surface temperature, 
changes the acoustical depth of the associated variation in 
Ti. This is illustrated in Fig. ^ for a selection of the models 
in Table 1; here Ti is plotted against the relative acousti- 
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Figure 6. The dashed line is Gas, the continuous line with error 
bars is G and the dot-dashed line G + d/uj 2 for the model of 
Table jl] indicated in each panel. To all the functions we have 
added a constant so that they match at the lowest frequency. 



T A 0 


Figure 7. Fi as function of relative acoustical depth r /tq for 
the ZAMS models in Table [l]. The continuous line is for Model 1, 
the dotted line for Model 3, the dashed line for Model 5 and the 
dot-dashed line for Model 7. 


els. We note that for all the models considered the filtered 
function G f agrees within the errors with the asymptotic 
G f as (which, for clarity, we do not show). As can be seen, for 
the models with M = 0.85 Mq, it is not possible to detect 
the signature of the Hell zone as given by G f even with fre¬ 
quency errors as small as 0.5 /xHz. However, for the rest of 
the models this frequency accuracy is sufficient. Also, from 
the figure it follows that G f (cu) has a larger amplitude for 
larger masses or later evolution stages, due to the increased 
strength of the He-II feature in IT (cf. Fig. |?j; see also Lopes 
et al. 1997). It can also be seen that in the fixed range of 
radial orders considered, G f contains fewer periods of the 
signal for M = 1.7 Mg than for the other masses. Indeed, it 
follows from equation (^), with A v ~ (2to) _1 , that in equa¬ 
tion © 2 lot ~ 27TOt/to. Thus for a fixed range in n, such 
as we have considered here, the signal arising from the sec¬ 
ond helium ionization zone generally oscillates less rapidly 
with decreasing relative depth t/tq. 


cal depth r/ro, where to is the acoustical radius of the star, 
obtained by integrating from r = 0 in equation ([Tt]). 

To determine the optimal p, we calculate the response 
function H(t) given by equation (^9|) for different values 
of p. Here the frequency interval and a as are computed for 
each model in Table |I| Then we choose p so as to keep the 
information from the second helium ionization zone, here 
defined at the point where Ti has a relative minimum. We 
have found that for Models 1 to 6 in Table [I] a value of p = 6 
is suitable. For Model 7 a value of p = 2 has been chosen, 
and for Model 8, p = 3. This decrease of p with increasing 
effective temperature clearly reflects the fact that the second 
helium ionization zone moves closer to the surface. 

Fig. ^ shows the functions G f ( u>) for the stellar mod- 
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3.2 Changes to the physics and parameters of the 
models 

We have considered several changes to the input physics and 
stellar parameters in order to analyse the sensitivity of the 
phase function to such quantities. Since, once the p-mode 
frequencies of a star are measured, the mean density can 
be accurately determined, it is interesting to compare mod¬ 
els with the same mean density, considering the effects of 
each modification separately. To do so, we have computed 
envelope models which have boundary conditions only at 
the surface, considering models with the same surface pa¬ 
rameters as those given in Table [l| For each of these eight 
models we have considered the modifications summarized 
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Figure 8. G f ( cu) for the model in Table h] indicated in each panel. 


in Table ||. Here a c is the mixing length parameter, X the 
hydrogen abundance, re the opacity, r c the radius at the 
base of the convection zone and EOS refers to the equation 
of state; all the models use the EFF formulation (Eggle- 
ton, Faulkner & Flannery 1973), except for the model corre¬ 
sponding to change 6 which uses the CEFF equation of state 
(see Christensen-Dalsgaard & Dappen 1992). For envelope 
models it is not possible to compute low-degree p-mode fre¬ 
quencies, and hence we cannot obtain G(u>) from a fit to 
equation (gj). However, as remarked previously, the differ¬ 
ences in the function G for the models corresponding to the 
small modifications considered in Table can be expected 
to be very similar to those in G as , which can be computed 
for envelope models. 

As an illustrative case, we shall show the results for 
just one case, an envelope model with the same global and 
surface parameters as Model 5 in Table In Fig.fi we show 


Table 2. Modifications to the envelope models. 



parameter 

magnitude 

notes 

1 

T e ft 

+100 K 


2 

Ot c 

-0.2 


3 

X 

-0.01 

at constant Z 

4 

logic K 

—0.2 in atmosphere 

at fixed r c 

5 

z 

+0.005 

at constant X 

6 

EOS 

CEFF 




v (mHz) 

Figure 9. The continuous line is Gas©) for an envelope model 
similar to Model 5 in Table []. The dashed line corresponds to the 
change in T e g, the dot-dashed line to the change in the mixing- 
length parameter a c , and the dotted line to the change in the at¬ 
mospheric opacity, for the changes listed in Table 2. The functions 
have been shifted by constants to match at a given frequency. 


Gas(w) for several modifications. The continuous line is for 
the reference model with errors corresponding to frequency 
errors of 0.5 fj, Hz. The dashed line corresponds to the change 
in T e ff, the dot-dashed line to the change in the mixing- 
length parameter q c , and the dotted line to the change in 
the atmospheric opacity. The functions G as for the changes 
in A', Z and the equation of state are not shown because 
they differ from G as for the reference model by less than 
the errors. This could be expected since these modifications 
do not change significantly the structure of the uppermost 
layers, contrary to those shown in the figure. Similar results 
are found for the other models in Table I. 

From Fig. 0, it follows that the function G computed 
from the observations can be used to constrain the structure 
of the uppermost layers. However, it is also clear that it is 
not possible to isolate one uncertainty from the rest, at least 
with an analysis of a single star. 

In Fig. [id, we show the differences <5Ga S (w) between the 
models with the changes indicated in Table 12 and the refer- 
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Figure 10. Differences in the filtered phase-function G'f between 
envelope models with the modifications indicated in Table ta and 
the reference model (corresponding to Model 5 in Table llT) for: 
change in T e g (continuous line), change in the mixing length pa¬ 
rameter a c (dotted line), change in the atmospheric opacities 
(long dashed line), change in X (dot-dashed line), change in Z 
(short-long dashed line) and change in the equation of state (short 
dashed line). 

ence model (an envelope model corresponding to Model 5 in 
Table [j]). As expected, <5G f is significant for the change in the 
equation of state or the envelope abundances, as a result of 
the corresponding changes in Fi in the second helium ioniza¬ 
tion zone; however, it is also important for modifications in 
the atmospheric opacities or the mixing-length parameter a c 
because these modifications change, for instance, the depth 
of the second helium ionization layer. This result is roughly 
similar to that found for the Sun (see Perez Hernandez & 
Christensen-Dalsgaard 1994b), although the relative impor¬ 
tance of each modification is different. 

For clarity we have not shown the errors in Fig. [Ifj but 
we note here that for frequency errors of 0.5 fi Hz, the mean 
error in G f is of 0.42 s and for 0.1 /rHz, the mean error in 
G f is 0.08 s. Hence, very accurate frequency determinations 
are required in order to obtain information from this func¬ 
tion, in particular concerning the envelope abundances or 
the equation of state. 


4 CONCLUSIONS 

By analysing low-degree (Z < 2) p-mode frequencies for the 
Sun we have shown that it is possible to fit an asymptotic 
expression which allows to separate the contribution of the 
upper layers as given by a function of frequency G(cu) from 
that of the interior. Moreover, it is possible to obtain a func¬ 
tion G f (w) that depends mainly on the layers around the 
second helium ionization zone. By applying the same tech¬ 
nique to p-mode frequencies of stellar models, we have found 


that the same separation is possible for main-sequence stars 
with masses between 0.85 Mq and 1.7 Mq. 

Of course, the main goal is to get information about 
stellar structure from these functions. When considering ob¬ 
served solar frequencies of low-degree modes, we have found 
that the function G(w) is determined with sufficient accu¬ 
racy to impose constraints on the structure of the uppermost 
layers of the Sun. Assuming frequency errors of 0.5 /rHz we 
found that the same can be true for other stars. Although 
different changes to the physics and parameters of the mod¬ 
els can lead to very similar G(u>), it might be possible to 
separate the effects of different uncertainties if several stars, 
belonging to the same cluster, are analysed; in this case it 
may be assumed, for example, that the stars share the same 
initial composition. 

In principle, analysis of the function G f could provide 
more significant constraints on the stellar models because 
it depends on the structure of the layers around the second 
helium ionization zone, where the physics is better under¬ 
stood than in the uppermost layers. However, very small er¬ 
rors are needed in order to obtain useful information from it. 
For instance, in the solar case, if only low-degree modes are 
used, frequency errors as small as those achieved by GOLF 
are required in order to determine the helium abundance 
in the solar envelope. Similar results are found for stars 
in the main sequence, for which errors substantially below 
0.5 /rHz are needed in order to impose significant constraints 
on the stellar models, for instance for the helium abundance 
or the equation of state. It is encouraging, therefore, that 
the COROT mission (Catala et al. 1995) aims at determin¬ 
ing frequencies with errors as small as 0.1 /.tHz. We note also 
that G f is more sensitive to modifications in the equilibrium 
structure for stars of mass greater than solar. 
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